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ABSTRACT 


It has been shown that rotational lines observed in the Horsehead nebula PDR 
are probably not caused by Z-C 3 H+, as was originally suggested. In the search 
for viable alternative candidate carriers, quartic force fields are employed here to 
provide highly accurate rotational constants, as well as fundamental vibrational 
frequencies, for another candidate carrier: 1 1 A' C 3 H”. The ab initio computed 
spectroscopic constants provided in this work are, compared to those necessary to 
define the observed lines, as accurate as the computed spectroscopic constants for 
many of the known interstellar anions. Additionally, the computed D e ff for CsH - 
is three times closer to the D deduced from the observed Horsehead nebula lines 
relative to Z-C 3 HA As a result, 1 1 A / CsH - is a more viable candidate for these 
observed rotational transitions and would be the seventh confirmed interstellar 
anion detected within the past decade and the first C n H~ molecular anion with 
an odd n. 

Keywords: astrochemistry — ISM: individual objects: Horsehead nebula — ISM: lines 
and bands — ISM: molecules — molecular data — radio lines: ISM 


1. Introduction 

Recent work by Huang et al. (2013) has questioned the attribution of lines observed 
in the Horsehead nebula photon-dominated- region (PDR) to Z-C 3 HA Quartic force fields 
(QFFs) computed from high-level ab initio quantum mechanical energies analyzed using 
perturbation theory at second order (Papousek & Aliev 1982) are known to produce 
highly accurate spectroscopic constants. Even though the B 0 computed by Huang et al. 
(2013) (11 262.68 MHz) is within 0.16% of the R-type rotational constant derived from the 
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observations by Pety et al. (2012) (11 244.9474 MHz), the computed D e of 4.248 kHz differs 
by 44.5% from the observed D value of 7.652 kHz. This “error” is more than an order of 
magnitude larger than any other error for a computed D e of a cation (using similar levels of 
theory) as compared to known high-resolution experimental data. Furthermore, the sextic 
distortion constant, H , differs by three orders of magnitude. As a result, it is unlikely that 
Z-C 3 H+ corresponds to the lines observed by Pety et al. (2012). 

This result motivates the question, “What is the carrier of these lines?” If these 
observed lines are, in fact, related to one another, certain inferences can be made about 
the molecular carrier. To match the rotational constants derived from the transition 
energies corresponding to the observed lines, the carrier is either linear or quasi-linear, 
almost certainly composed of three carbon atoms as well as a single hydrogen atom, and 
closed-shell since there are no splittings in the lines as required for the rotational spectra 
of open-shell molecules (McCarthy 2013). All of these criteria are, in fact, met by l- C 3 H + , 
but this cation’s difference between observational and high-accuracy theoretical rotational 
constants, especially the D constant, discussed above and by Huang et al. (2013), probably 
rules it out. As a result, the quasi-linear anion, 1 1 A / /-C 3 H~, remains as the most likely 
candidate carrier of the Horsehead nebula PDR rotational lines of interest. 

Even though the most stable singlet isomer of C 3 H~ is the cyclic form, c-C 3 H~, the 
barrier to isomerization is high enough ( 45 kcal mol -1 ) for the quasi-linear C s isomer to be 
kinetically stable (Lakin et al. 2001). Various mechanisms for interstellar synthesis of this 
anion are possible (Larsson et al. 2012; Senent & Hochlaf 2013) and are probably related 
to those responsible for the creation of the related C 2 „H _ for n = 2 — 4 anions previously 
detected in the interstellar medium (ISM) (McCarthy et al. 2006; Cernicharo et al. 2007; 
Briinken et al. 2007b). Furthermore, radical C 3 H in both the linear and cyclic forms has 
also been detected in the ISM (Thaddeus et al. 1985; Yamamoto et al. 1987) suggesting the 
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possible interstellar existence of the anion. 


Additionally, CsH“ is of astronomical interest since it has been computationally shown 
by Fortenberry (2013) to possess not only a rare dipole-bound singlet excited electronic state 
(the 2 l A' state) but also an even more rare valence excited state (1 1 A") below the electron 
binding or electron detachment energy. In fact, the valence electronically excited state is 
the only such state thus far proposed to exist for an anion of this size which also contains 
only first-row atoms (Fortenberry & Crawford 201 lb, a; Fortenberry 2013). Furthermore, 
anions have been proposed as carriers of some diffuse interstellar bands (DIBs) (Sarre 2000; 
Cordiner & Sarre 2007; Fortenberry et al. 2013a), and C3H - has unusual properties that 
may be of importance to the DIBs. 


2. Computational Details 

The spectroscopic constants and fundamental vibrational frequencies of 1 l A' /-C3H" 
are computed through the established means of QFFs (Huang & Lee 2008). Starting from a 
restricted Hartree-Fock (RHF) (Schemer et al. 1987) coupled cluster (Lee & Scuseria 1995; 
Shavitt & Bartlett 2009; Crawford & Schaefer 2000) singles, doubles, and perturbative 
triples [CCSD(T)] (Raghavachari et al. 1989) aug-cc-pV5Z (Dunning 1989; Kendall et al. 
1992; Dunning et al. 2001) geometry further corrected for core correlation effects from the 
Martin-Taylor (MT) basis set (Martin & Taylor 1994), a grid of 743 symmetry-unique 
points is generated. Simple-internal coordinates for the bond lengths and Z H— C— C are 
coupled to linear LINX and LINY (Allen & coworkers 2005) coordinates exactly as those 
defined in Fortenberry et al. (2012b) for HOCO + . Displacements of 0.005 A for the bond 
lengths, 0.005 rad for the bond angle, and 0.005 for the LINX and LINY coordinates and 
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the associated energies computed at each point define the QFF, which is of the form: 

f = 2 ^ ] FijAiAj + — N ] -FjfcjAjAjAfc + — ^ ^ FikjiAiA.jA.kAi, (1) 

ij ijk ijkl 

where Aj are the displacements and F l3 . are force constants (Huang & Lee 2008). 

At each point, CCSD(T)/aug-cc-pVXZ (where X = T,Q, 5) energies are computed and 
extrapolated to the complete basis set (CBS) limit via a three-point formula (Martin & Lee 
1996). Additionally, energy corrections are made to the CBS energy for core correlation 
and for scalar relativistic effects (Douglas & Kroll 1974). The resulting QFF is denoted 
as the CcCR QFF for the CBS energy, core correlation correction, and scalar relativistic 
correction, respectively, (Fortenberry et al. 2011). The augmented Dunning basis sets have 
been shown by Skurski et al. (2000) to be reliable for computations of anionic properties. 
An initial least-squares-fit of the CcCR energy points leads to a minor transformation of 
the reference geometry such that the gradients are identically zero. This geometry and 
the resulting force constants are then employed in the rovibrational computations. All 
electronic structure computations make use of the MOLPRO 2010.1 quantum chemical 
package (Werner et al. 2010), and all employ the Born-Oppenhcimer approximation making 
the QFFs identical for the isotopologues. 

The QFF is fit from the 805 redundant total energy points with a sum of squared 
residuals on the order of 3 x 10~ 17 a.u. 2 Cartesian derivatives are then computed from the 
QFF with the INTDER program (Allen & coworkers 2005). From these, the SPECTRO 
program (Gaw et al. 1991) employs second-order vibrational perturbation theory (VPT2) to 
generate the spectroscopic constants (Papousek & Aliev 1982) and vibrational frequencies 
(Mills 1972; Watson 1977). After transforming the force constants into the Morse-cosine 
coordinate system so that the potential possesses proper limiting behavior (Dateo et al. 
1994; Fortenberry et al. 2013b), vibrational configuration interaction (VCI) computations 
with the MULTIMODE program (Carter et al. 1998; Bowman et al. 2003) also produce 


5 


vibrational frequencies. The VCI computations make use of the same basis set and 
wavefunction construction as that utilized by Fortenberry et al. (2012b, a) in similar 
quasi-linear tetra-atomic systems. 


3. Discussion 

The force constants computed in this study are listed in Table 1. The CcCR 
geometrical parameters and spectroscopic constants are given in Table 2 for both 1 1 A' 
l- C3H - and the deuterated isotopologue. The equilibrium dipole moment is computed 
with respect to the center-of-mass with CCSD(T)/aug-cc-pV5Z to be 2.16 D. The C— C— C 
R a vibrationally-averaged bond angle is nearly collinear at 174.540° while the H— C— C 
vibrationally- averaged bond angle is 109.491°. These values are in line with those computed 
by Lakin et al. (2001). As has been discussed by Fortenberry (2013) for C 3 H~, the Ci 
carbon atom adjacent to the hydrogen atom shown in Figure 1 is a carbene-type carbon 
containing a lone pair which leads to a longer Ci— C2 bond length compared to the shorter 
C2 — C3 bond length. Even though this result differs from the CCSD(T) results from Lakin 
et al. (2001), their reported CASSCF and HF results give bond lengths similar to ours 
leading us to conclude that the CCSD(T) C— C bond lengths are mislabeled in the paper 
by Lakin et al. (2001). The vibrationally-averaged geometrical parameters change slightly 
upon deuteration. Similar bond angles of the heavy atoms have been computed for the 
frans-HOCO + , HOCS + , and HSCO + systems (Fortenberry et al. 2012b, a) with very good 
agreement present for known experimental data. 

The most notable values in Table 2 are the rotational constants and the quartic 
centrifugal distortion (D-type) constants. For 1 1 A' Z-C 3 H~, the Bo rotational constant is 
11 339.66 MHz while Co is 11 087.35 MHz. The equilibrium constants are slightly larger, 
but both sets are in reasonable agreement with those computed by Lakin et al. (2001). The 
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.D-type constants have not been vibrationally averaged, and Dj , most prominently, is 4.954 
kHz. 


Direct comparison between these values and those deduced from the Horsehead nebula 
PDR spectrum observed by Pety et al. (2012) is not possible since the isomer of C 3 H“ of 
interest here is not perfectly linear. Pety et al. (2012) assume a linear structure in order to 
fit the effective rotational constant, B e ff, and the effective centrifugal distrotion constant, 
D e fj and use the second-order fitting equation, 

vj+i^j = 2 B(J + 1) — 4 D[J + l) 3 , (2) 


to compute the affiliated rotational constants. C 3 H is non-linear and requires the following 
related equation from McCarthy et al. (1997): 


U/+i->j — (B + C)(J + 1) — < 4Dj + 


( B ~ Cf 


A - 


(B+C) 


(J + If, 


( 3 ) 


with the assumption that K = 0 forcing c = 8. As such, we can set Equation 2 equal to 
Equation 3. The (J + 1) term in Equation 3 is equal to 2 -B e //, and the (J + l) 3 term in 
Equation 3 is equal to 4 D e ff. Using the CcCR computed A 0 , B 0 , C 0 , and Dj values, where 
Dj is the only equilibrium constant, B e ff is computed to be 11 213.51 MHz, and D e ff is 
8.795 kHz. Hence, direct comparison between the CcCR C 3 H~ rotational constants and 
those obtained from the lines observed by Pety et al. (2012) is possible. 


The second-order fit of the lines observed by Pety et al. (2012) indicates that the carrier 
must have a B/B e ff that is very close to 11 244.9474 MHz and a D/D e ff that is 7.652 kHz. 
The B e ff computed with the A 0 , H 0 , and C 0 rotational constants by the above approach 
is very close, off by 31.44 MHz or 0.28%. This is roughly the same difference between the 
observed B and that of Z-C 3 H + (Huang et al. 2013) and is typical of the accuracies that 
can be expected when comparing QFF computations with the CCSD(T) level of theory to 
experiment as enumerated in Table 3. However, the 8.795 kHz D e ff for 1 l A' Z-C 3 H” is 
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much closer to the 7.652 kHz D derived from the lines observed by Pety et al. (2012) in the 
Horsehead nebula than the linear cation (Huang et al. 2013). Even so, this D e ff of 8.795 
kHz differs from the observation by 1.14 kHz or 14.93%. Replacing Ho, Bq, and Cq with A e , 
B e , and C e in Equation 3 lowers the C 3 H - -D e // to 8.366 kHz, a difference of 0.714 kHz or 
9.34% from that determined by Pety et al. (2012). 

Table 3 provides some insight into the accuracies that can be expected for calculated 
rotational constants of similar molecules. The quasi-linear cations listed below CsH” show 
strong correlation between the computed B e ff from the use of B 0 and C 0 and the B e ff 
derived from the various experiments. Additionally, the D e ff values computed the same 
way with the equilibrium Dj also show good, albeit not as strong, correlation between 
theory and experiment. Unfortunately, CaH” has errors that are larger than this. However, 
this probably results from a combination of basis set incompleteness and higher-order 
correlation effects. Even though aug-cc-pVXZ basis sets used at the CCSD(T) level of 
theory have been shown to be effective in the computation of anionic properties (Skurski 
et al. 2000; Fortenberry & Crawford 2011b), higher-order properties such as the ZUtype 
constants are more susceptible to even the smallest errors. This is clear in the cations, as 
well, where the D e ff values are not as accurate as the -B e // values. 

The known interstellar anions and the related C 2 H~ system, which has not yet been 
detected in the ISM, are linear and have B and D computed directly, most often as B e 
and D e . Note that the theoretical rotational constants are not as accurate for the anions 
as they are for the cations. Most notably, the B e /Bo and D e values computed with a 
CCSD(T)/aug-cc-pCV5Z cubic force field for C 5 N - by Botschwina & Oswald (2008) are 
directly used in the identification of this anion in the ISM (Cernicharo et al. 2008). As 
listed in Table 3, agreement between computed B values and that necessary to match 
the observed rotational lines actually worsens when Bq is used instead of B e , more than 


doubling the percent error. This is the same behavior as what is currently found for C 3 H - . 
Additionally, the D e percent error for C 5 N - , as compared to observation, is 9.1%, almost 
exactly what it is for C 3 H - when using the equilibrium rotational constants. The force 
held employed by Botschwina & Oswald (2008) also includes core correlation like the CcCR 
QFF. Hence, the present rotational constants are in the same accuracy range for CsH - as 
those used to detect CsN~ in the ISM. 

Furthermore, Senent & Hochlaf (2010) compute B e for C 4 H~ to be in error compared 
to experiment by nearly twice as much as B e ff for CsH - utilizing A 0l B 0l and Co in the 
formulation here. The C 2 H" B e is in error compared to experiment by roughly an equivalent 
amount as B e ff for C 3 1 1 . Also, 1), for Cfd R and C'sII is in the same error range as the 
present C 3 H~ -D e //- C 3 N _ behaves similarly. The CCSD(T)/aug-cc-pCV5Z B e and D e for 
CN - are much more accurate than any other system shown here, but it has fewer electrons 
and fewer degrees of freedom than any other. Regardless, our accuracies in the computation 
of B e ff and -D e // for 1 1 A / C 3 H - are comparable to those obtained for other interstellar 
anions. Hence, CsH~ is, perhaps, the most viable candidate for the carrier of the rotational 
lines observed in the Horsehead nebula previously attributed to Z-C 3 H+ (Pety et al. 2012). 

Comparison of the sextic distortion constant, H e ff, is not as straightforward. There is 
a dearth of data on how the computation of this value for anions compares to experiment. 
Hj , which is an equilibrium value, is not exactly H e ff , but they are probably related. Even 
though H obtained by Pety et al. (2012) is 560 mHz and Hj for C 3 H~ is 3.344 mHz, this is 
an order of magnitude closer agreement than this same H compared to the H e for /-C 3 H + , 
0.375 mHz (Huang et al. 2013). Additionally, the same basis set and correlation errors for 
anions that affect the calculation of D will be present for H. As a result, we can only say 
here that as far as H is concerned for comparison to the lines observed in the Horsehead 
nebula by Pety et al. (2012), 1 1 A / CsH~ is a better candidate than Z-C 3 HA 
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The harmonic and anharmonic vibrational frequencies for both 1 1 A' /-C 3 II and 
/-CaD - are given in Table 4. Positive anharmonicities are present in both isotopologues 
for the 1/5 Ci— C 2 — C 3 bending and the z / 6 torsional modes. VPT2 and VCI produce 
fundamental vibrational frequencies from the CcCR QFF that are quite consistent. The 
largest deviation between the methods, 1.0 cm -1 , is found for the z / 4 H— Ci— C 2 bending 
mode. Comparison of the C 3 I R CcCR QFF vibrational frequencies, whether computed 
using VPT2 or VCI, to those computed by Lakin et al. (2001) is roughly consistent for 
iq-iq. The u 5 anharmonic frequencies differ by more than 50 cm -1 , though the CJ 5 harmonic 
frequencies are very similar (i.e., the difference in the 1/5 fundamental frequency is mostly 
due to differences in the anharmonic correction). The torsional mode is nearly identical 
between the two studies, though in this case the harmonic frequencies differ by more than 
50 cm -1 . It is hoped that the present QFF computations of the fundamental vibrational 
frequencies provided here will assist in the characterization of this anion in current and 
future studies of the ISM or simulated laboratory experiments at infrared wavelengths in 
addition to studies in the sub-millimeter spectral region. 


4. Conclusions 

Since the link between /-C 3 II 1 and the lines observed in the Horsehead nebula PDR 
by Pety et al. (2012) has recently been strongly questioned by Huang et ah (2013), another 
viable candidate is necessary. The rotational lines seem to require a closed-shell quasi-linear 
structure composed of three carbon atoms along with a hydrogen atom. 1 l A' CsH~ appears 
to be the most likely candidate. Here, the CcCR QFF has determined a B e ff for this 
anion to be in error by 0.28% from that required to fit the observed lines. The use of the 
equilibrium rotational constants fortuitously lowers the error to 0.11%. However, the error 
reduction and error magnitudes themselves are in line with the computed CsN - rotational 


10 


constants used in its interstellar detection. Additionally, the discrepancy between the A e , 
B e , and C e computed C 3 H~ D e ff and the As// deduced from the observed interstellar 
rotational lines is similar to the D e errors for C 4 H - , C 3 N _ , and CsN - and less than that 
of CgH - , which are all reported for CCSD(T) computations, i.e. similar levels of theory. 
Hence, the consistency of the errors for C 3 H - with other anions previously detected in the 
ISM coupled with its matching the required spectral criteria make this anion the strongest 
candidate carrier for the Horsehead nebula PDR lines and, potentially, the most recent 
anion detected in the ISM. It would also be the first detected interstellar odd-numbered 
monohydrogen carbon chain anion. 
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Fig. 1. — CcCR equilibrium geometry of 1 l A' Z-C 3 H . 
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Table 1: The simple- internal CcCR QFF Quadratic, Cubic, and Quartic Force Constants (in 
mdyn / A n -rad m ) 0 for Z-C 3 H~. 


Fii 
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-F 431 

0.0711 

Thin 

318.24 

7*4432 

0.24 

7*5531 

0.12 
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7*4433 
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7*5553 

0.12 

F 5 1 

0.070 974 

F531 

-0.0597 

F3321 

-0.80 

7*5311 

0.09 

775554 

- 0.12 

F52 

0.081 130 

7*532 

-0.2190 

F3322 

-0.54 

7*5321 

0.21 

7*5555 

1.95 

F53 

0.069 481 

7*533 

-0.0030 

F3331 

-0.82 

7*5322 

0.54 

7*6611 

0.14 

7*54 

0.064 059 

7*541 

-0.0774 

7*3332 

0.46 

77)331 

-0.04 

7*6621 

-0.49 

7*55 

0.404 485 

7*542 

0.0376 

7*3333 

145.05 

77)332 

-0.07 

7*6622 

0.64 

7*66 

0.168 044 

F543 

-0.0694 

7*4111 

-0.17 

7*5333 

-0.46 

7*6631 

0.03 

Fin 

-64.7214 

7*544 

-0.1425 

F4211 

- 0.02 

7*5411 

0.14 

7*6632 

0.06 

Fn\ 

0.5759 

F551 

-0.4284 

F4221 

-0.27 

7*5421 

0.12 

7*6633 

-0.14 

F221 

-3.2972 

7*552 

-0.9210 

F4222 

-0.32 

7*5422 

-0.49 

7*6641 

-0.07 

F222 

-43.3783 

7*553 

-0.0900 

771311 

0.15 

77)431 

-0.05 

7*6642 

-0.04 

F311 

0.0659 

-F554 

0.0071 

7*4321 

-0.03 

77)432 

0.31 

773643 

- 0.10 

F321 

-0.3042 

F555 

-0.1839 

7*4322 

-0.46 

7*5433 

-0.09 

773644 

0.08 

F322 

-0.0840 

Fmi 

-0.1710 

F4331 

-0.06 

7*5441 

0.01 

7*6651 

-0.04 

F33 1 

0.1287 

Fffi2 

-0.3467 

F4332 

-0.31 

77)442 

0.14 

7*6652 

0.02 

7*332 

0.2601 

7*663 

-0.0476 

F 4333 

-1.42 

77)443 

0.24 

7*6653 

-0.03 

7*333 

-28.7819 

7*664 

0.0133 

7*4411 

-0.63 

77)444 

-0.06 

7*6654 

- 0.10 

7*411 

-0.2017 

7*665 

-0.0708 

7*4421 

1.45 

7*5511 

0.57 

7*6655 

0.23 

F42I 

0.3282 



7*4422 

-1.71 

7*5521 

0.14 

7*6666 

0.86 

F422 

-0.6518 



F4431 

-0.17 

7*5522 

1.56 




a l mdyn =10 8 N; n and m are exponents corresponding to the number of units from the type of 

modes present in the specific force constant. 
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Tabic 2: The Zero-Point (R a vibrationally- averaged) and Equilibrium Structures, Rotational 
Constants, CCSD(T)/aug-cc-pV5Z Dipole Moment, Vibration- Rotation Interaction Con- 
stants, and Quartic and Sextic Distortion Constants of 1 1 R / Z-C 3 H - and the deuterated 
form with the CcCR QFF. 



c 3 h- 

Previous 42, 

c 3 d- 

r 0 (Cl— H) 

1.119 438 A 


1.116 446 A 

ro(Ci-Ca) 

1.351 595 A 


1.351 753 A 

ro(C 2 -C 3 ) 

1.282 845 A 


1.282 620 A 

Z 0 (H-C 1 — C 2 ) 

109.491° 


109.530° 

^o(Ci-C 2 -C 3 ) 

174.540° 


174.643° 

A 0 

529 134.2 MHz 


295 539.6 MHz 

B o 

11 339.66 MHz 


10 626.03 MHz 

Co 

11 087.35 MHz 


10 238.74 MHz 

Bj 

4.954 kHz 


4.544 kHz 

Djk 

0.702 MHz 


0.316 MHz 

Dk 

217.543 MHz 


94.897 MHz 

d\ 

-0.112 kHz 


-0.253 kHz 

d 2 

-0.023 kHz 


-0.052 kHz 

Hj 

3.344 mHz 


16.516 mHz 

Hjk 

3.221 Hz 


2.151 Hz 

HkJ 

-3.229 kHz 


-0.745 kHz 

Hk 

358.867 kHz 


90.731 kHz 

«1 

0.132 mHz 


0.634 mHz 

H 2 

0.203 mHz 


0.612 mHz 

h 3 

0.037 mHz 


0.133 mHz 

Tn.n.nn 

-873.001 MHz 


-380.872 MHz 

r bbbb 

-0.021 MHz 


-0.021 MHz 

Tcccc 

-0.019 MHz 


-0.017 MHz 


-2.766 MHz 


-1.619 MHz 


-0.081 MHz 


0.319 MHz 


-0.020 MHz 


-0.018 MHz 

$ aaa 

355 640.661 Hz 


89 988.504 Hz 

®bbb 

0.001 Hz 


0.004 Hz 

4*ccc 

0.000 Hz 


0.001 Hz 

®aab 

390.158 Hz 


703.204 Hz 

®a.bb 

4.265 Hz 


3.112 Hz 

aac 

-3 614.354 Hz 


-1 445.590 Hz 

®bbc 

0.000 Hz 


0.001 Hz 

acc 

-0.271 Hz 


0.151 Hz 


0.001 Hz 


0.002 Hz 

*obc 

4.570 Hz 


3.618 Hz 

cc A 1 

27 922.5 MHz 


11 662.9 MHz 

2 

-725.5 MHz 


-917.5 MHz 

3 

484.8 MHz 


170.2 MHz 

a A 4 

-35 092.1 MHz 


-16 226.3 MHz 

c* A 5 

-3 103.1 MHz 


-4 597.5 MHz 

a A 6 

12 333.4 MHz 


9 042.1 MHz 

a B 1 

4.2 MHz 


6.9 MHz 

a B 2 

83.5 MHz 


77.2 MHz 

c B 3 

45.1 MHz 


40.3 MHz 

ol B 4 

-12.0 MHz 


-8.4 MHz 

a B 5 

-47.1 MHz 


-48.4 MHz 

ol B 6 

-48.6 MHz 


-45.9 MHz 

a C 1 

14.8 MHz 


18.3 MHz 

OL C 2 

78.6 MHz 


70.1 MHz 

c c 3 

38.4 MHz 


39.3 MHz 

a C 4 

16.0 MHz 


12.9 MHz 

a C 5 

-16.1 MHz 


-15.2 MHz 

OtP 6 

-78.5 MHz 


-69.8 MHz 

r e (Ci-H)'> 

1.106 939 A 

l.no A 

- 

r e(Ci C 2 ) 

1.349 832 A 

1.289 A 

- 

r e (C 2 C 3 ) 

1.281 900 A 

1.363 A 

- 

Z e (H-Ci-C 2 ) 

109.529° 

109.2° 

- 

Z e (C 2 -C 3 -C 4 ) 

174.571° 

171.2° 

- 

A e 

530 044.3 MHz 

524.5 MHz 

295 106.5 MHz 

Be 

11 352.05 MHz 

11.2 MHz 

10 636.73 MHz 

Ce 

11 114.02 MHz 

10.9 MHz 

10 266.68 MHz 


2.16 D 

- 

- 

Mcc 

1.63 D 

- 

- 

Vy 

1.41 D 

- 

- 


a CCSD(T)/aug-cc-pVQZ QFF results from Lakin et al. (2001). 

^The equilibrium geometries are identical among isotopologues from the use of the Born-Oppenheimer approximation. 
c The C 3 H — coordinates (in A with the center-of-mass at the origin) used to generate Born-Oppenheimer dipole moment components are: H, 
1.733414, -0.910473, 0.000000; Ci, 1.276456, 0.098036, 0.000000; C 2 , -0.069613, -0.016965, 0.000000; C 3 , -1.352424, -0.004605, 0.000000. 
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Tabic 3: Errors in the computation of B e ff (in MHz) and D e ff (in kHz) for various small 
molecules. 


Molecule 

Experiment 

B / B eff 

Theory 

% Error 

D / D eff 

Experiment Theory 

% Error 

c 3 h-“ 

11244.9474 

11213.51 

0.28% 

7.652 

8.795 

14.93% 

C 3 H~ b 

11244.9474 

11233.04 

0.11% 

7.652 

8.366 

9.34% 

HSCO +c 

5636.866 

5637.60 

0.01% 

3.1 

3.116 

0.53% 

HOCO +d 

10691.58265 

10705.44 

0.13% 

4.580576 

4.511 

1.52% 

HOCS +e 

5726.66011 

5730.22 

0.06% 

1.064 

1.107 

4.00% 

c 2 h-' 

41639.20 

41781.0 

0.34% 

0.09697 

0.0946 

2.44% 

c 2 h-^ 

41639.20 

41614.0 

0.06% 




C 4 H-9 

4654.9449 

4653.9 

0.02% 

0.5875 

0.55 

6.4% 

CiR- h 

4654.9449 

4625.6546 

0.63% 




c 6 h-* 

1376.86298 

1376.9 

0.00% 

0.03235 

0.0270 

16.5% 

c 8 h-* 

583.30404 

583.2 

0.02% 

0.0042 

0.0033 

17% 

CN~ fc 

56132.7562 

56152 

0.03% 

186.427 

185 

0.77% 

C 3 N _i 

4851.62183 

4850 

0.03% 

0.68592 

0.628 

8.44% 

C 5 N- m 

1388.860 

1387.8 

0.08% 

0.033 

0.0300 

9.1% 

C 5 N~ n 

1388.860 

1386.2 

0.19% 





“This work with vibrationally averaged data and the observational lines described by Pety et al. 

( 2012 ). 

^This work from the equilibirum rotational constants and the observational lines described by 
Pety et al. (2012). 

c CcCR QFF data computed from vibrationally averaged A, B, and C constants given in 
Fortenberry et al. (2012a) and experimental data from Ohshima & Endo (1996). 
rf CcCR QFF vibrationally averaged data from Fortenberry et al. (2012b) and experimental data 
from Bogey et al. (1988). 

e CcCR QFF vibrationally averaged data from Fortenberry et al. (2012a) and experimental data 
from McCarthy Sz Thaddeus (2007). 

■^CcCR QFF equilibrium values (first line) and Bq (second line) from Huang & Lee (2009) and 
experimental data from Briinken et al. (2007a). 

9 CCSD(T)/cc-pVTZ B e and D e as well as experimental data from Gupta et al. (2007). 
/l RCCSD(T)/aug-cc-pVQZ B e from Senent & Hochlaf (2010) with the experimental B e fj again 
from Gupta et al. (2007). 

*CCSD(T)/cc-pVTZ B e , CCSD(T)/cc-pVDZ D e , and experimental data from McCarthy et al. 
(2006). 

■ ? CCSD(T)/cc-pVTZ B e , SCF/DZP D e , and experimental data from Gupta et al. (2007). 
fc CCSD(T)/aug-cc-pCV5Z B e , CCSD(T)/aug-cc-pCVQZ D e , and experimental data from 
Gottlieb et al. (2007). 

z CCSD(T)/aug-cc-pCV5Z B e and D e from Kolos et al. (2008) (ABq is reported to be 0.606 MHz 
giving a % error of about 0.02%) and experimental data from Thaddeus et al. (2008). 
m CCSD(T)/aug-cc-pCV5Z B e and CCSD(T)/aug-cc-pVQZ D e from Botschwina & Oswald 
(2008) with experimental data from Cernicharo et al. (2008). 

n CCSD(T)/aug-cc-pCV5Z Bo from Botschwina & Oswald (2008) with experimental data from 
Cernicharo et al. (2008). 19 


Tabic 4: The C 3 H and C 3 D CcCR QFF harmonic, VCI, and VPT2 fundamental vibra- 


tional frequencies in cm 


Mode 

Description 

Harmonic 

c 3 h- 

VCI 

VPT2 

Previous 0, 

Harmonic 

C 3 H- 

Anharm. 

Harmonic 

c 3 d~ 

VCI 

VPT2 

v\ (a') 

Ci— H stretch 

2881.9 

2714.4 

2713.9 

2863 

2723 

2122.9 

2036.4 

2035.5 

^2 (a') 

C 2 — C 3 stretch 

1843.9 

1804.3 

1804.4 

1831 

1828 

1832.9 

1796.5 

1796.5 

^3 (o') 

Ci— C 2 stretch 

1117.1 

1108.0 

1107.9 

1091 

1120 

1112.0 

1100.9 

1101.0 

j/ 4 (a') 

H-C 1 -C 2 bend 

1037.8 

1012.1 

1011.1 

1002 

1022 

817.0 

803.8 

802.7 

"s(a') 

C 1 -C 2 -C 3 bend 

406.7 

419.4 

418.9 

393 

368 

379.1 

382.4 

381.9 

u 6 (a") 

torsion 

281.0 

296.8 

296.1 

349 

297 

278.9 

286.7 

286.1 


a CCSD(T)/aug-cc-pVQZ QFF results from Lakin et al. (2001). 
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